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In a recent paper, we analyzed the properties of a new kind of spherical wavelets (called needlets) 
for statistical inference procedures on spherical random fields; the investigation was mainly 
motivated by applications to cosmological data. In the present work, we exploit the asymptotic 
uncorrelation of random needlet coefficients at fixed angular distances to construct subsampling 
statistics evaluated on Voronoi cells on the sphere. We illustrate how such statistics can be 
used for isotropy tests and for bootstrap estimation of nuisance parameters, even when a single 
realization of the spherical random field is observed. The asymptotic theory is developed in 
detail in the high resolution sense. 
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1. Introduction 

Strong empirical motivations have prompted rising activity in the statistical analysis of 
spherical random fields over the last few years, mostly in connection with the analy- 
sis of the cosmic microwave background (CMB) radiation. The first CMB maps were 
provided by the NASA mission COBE in 1993 and led to the Nobel prize for physics 
for J. Mather and G. Smoot in 2006. Much more refined observations were provided 
by the satellite mission WMAP in 2003/2006. Further improvements are expected in 
the next few years, in view of the forthcoming launch of the European Space Agency 
mission Planck. Cosmic microwave background data have been collected by many other 
remarkable balloon-borne experiments; a list with links to some data sets can be found 
on http : / / www. fisica. uniroma2 . it / ~ cosmo /. 

The analysis of CMB data provides a sort of goldmine of new challenges for statistical 
methodology; see, for instance, Dodelson [8] for a review. Many of these challenges can 
be addressed by some forms of spherical wavelets; to list just a few, we recall testing 
for non-Gaussianity (Cabella et al. [6], Vielva et al. [29], McEwen et al. [17], Jin et al. 
[13]), component separation (Mouddcn et al. [18]), foreground subtraction (Hansen et al. 
[11]), point sources detection (Sanz et al. [25]), cross-correlation with large scale structure 
data (Pietrobon et al. [22], McEwen et al. [16]), searching for asymmetries/directional 
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features in CMB (Wiaux et al. [32], Vielva et al. [30]) and many others. The importance 
of wavelets in this environment is easily understood: on one hand, nearly all predictions 
from theoretical physics for the behaviour of the CMB field are presented in Fourier space, 
where the orthogonality properties make many difficult problems more tractable. On the 
other hand, spherical CMB maps are usually observed with gaps due to the presence 
of foreground radiation such as the emissions by the Milky Way and other galaxies. 
The double localization properties of wavelets thus make them a very valuable asset for 
CMB data analysis. Because of this, indeed very many physical papers have attempted 
to use wavelets in cosmology for a variety of different problems. Many of these papers 
provide insightful suggestions and important experimental results; however, the focus has 
typically been on physical data and few have investigated the statistical foundations of the 
proposed procedures. In most cases, rather than proposing anew procedures suitable for 
spherical random fields, it has been very common to rely on tangent plane approximations 
to adapt to the sphere standard wavelet constructions on the plane (an exception is 
provided by Sanz et al. [25]; see also Antoine and Vandergheynst [1], Antoine et al. [2] 
and Wiaux et al. [31] for a nice group-theoretical construction). 

Nccdlets are a new kind of second-generation spherical wavelets, which were introduced 
into functional analysis by Narcowich, Petrushev and Ward [19, 20]; they can be shown 
to make up a tight frame with excellent localization properties in both the real and the 
harmonic domains. In Baldi et al. [4], their properties for the analysis of random fields 
were established; in particular, it was shown that a major feature of random nccdlct 
coefficients is their asymptotic uncorrelation at large frequencies j for any fixed angular 
distance (see Baldi et al. [3] for an analogous result on the circle). Of course, in the 
Gaussian and isotropic cases, this property implies that the random spherical needlets 
behave asymptotically as an i.i.d. array. This surprising result clearly opens the way 
to the implementation of several statistical procedures with an asymptotic justification. 
The meaning of asymptotics in this framework should be understood with great care. 
It should be stressed that we are considering a single observation (our universe) of a 
mean square continuous and isotropic random field on a fixed domain. Our asymptotic 
theory is then entertained in the high-resolution sense (compare Marinucci [15]), that is, 
as higher and higher frequency data become available for statistical analysis. 

In this paper, we build on this essential feature of the random needlet coefficients to 
propose new statistical procedures and to provide bootstrap estimates for the asymptotic 
variance of existing techniques. The main idea of our work can be described as follows: We 
partition the sphere S 2 into disjoint subsets with roughly the same dimension, in a sense 
to be made rigorous later. It is then possible to evaluate nonlinear statistics on each of 
these subregions separately. These subregions will constitute a tessellation of the spherical 
surface composed by Voronoi cells associated with a suitable e-net. The geometry of 
these Voronoi cells on the sphere is interesting by itself and plays an important role in 
our results. 

We show below that, in the high-resolution sense, these statistics converge asymptot- 
ically to a sequence of independent realizations with the same law as the corresponding 
ones evaluated on the full sphere. It is then immediate to exploit this result to obtain 
computationally feasible approximations of sample variance for nonlinear functionals (in 
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Baldi et al. [4] estimation of the normalizing factors was left open). It is also quite 
straightforward to exploit our construction to derive tests for statistical isotropy on the 
sphere. The latter issue is indeed of great importance for cosmological data analysis: in- 
deed the single most surprising result from the first releases of WMAP data in 2003/2006 
was the apparent presence of an asymmetry in CMB radiation (see, for instance, Hansen 
et al. [11])- These findings have sparkled an impressive amount of empirical research 
over the last few years, but the results are still inconclusive, partly due to the lack of 
widely accepted statistical procedures to tackle this issue. The actual discovery of an 
asymmetry in cosmological radiation might call for revolutionary advances in theoretical 
physics, possible explanations ranging from the higher-dimension non-trivial topological 
structure of the observed universe to rotating solutions of Einstein field equations. 

To some extent, our approach can be viewed as a simple form of subsampling in the 
sense of Politis et al. [23]. To the best of our knowledge, subsampling techniques have 
so far been considered mainly in the presence of the usual large sample asymptotics, 
the main instrument to establish asymptotic properties being some mixing properties 
as the observations increase. In the present circumstances, no mixing properties can be 
advocated, as we assume we are observing a single realization of a random field on a 
bounded and compact domain. Once again, then, we consider our results to be a quite 
surprising consequence of the peculiar properties of the needlets construction. 

The structure of the paper is as follows: in Section 2 we review some basic features 
of the needlets and their properties in the analysis of isotropic random fields. Section 
3 discusses the geometry of the sphere and, more precisely, the separation properties 
of Voronoi cells and cubature points inside them. In Section 4 we provide our main 
theoretical results, which are made possible by careful correlation inequalities that are 
justified in the Appendix. Section 5 discusses the results, their statistical applications 
and possible routes for further research. 

2. A review on needlets 

In this section, we review very briefly a few basic features of needlet construction. For 
more details, we refer to Narcowich, Petrushev and Ward [19] and Baldi et al. [4]. Let us 
denote by {Kz m } m =-z....,z the spherical harmonics (see [28]), that form an orthonormal 
basis for L2(S 2 ,/i). The following decomposition holds: 

L 2 (§ 2 , M ) = 0ffz, 

1=0 

where the H^s are the finite dimensional spaces of L 2 spanned by the Ith spherical 
harmonics and [i is Lebesgue measure on the sphere. We define also the space of the 
restrictions to S 2 of the polynomials of a degree not greater than I as 

/ 

m=0 



Subsampling needlet on the sphere 441 

Now let <j> be a C°° non-increasing function supported in |£| < 1, such that 1 > <f>(£) > 
and 0(£) = 1 if |£| < B' 1 , and define: 

so that 

vki>i, E &2 (ij) = 1 ' B>1 - C 1 ) 

i 

Define the projection operator 

L l((x,y)) = E Y lm {x)Yi m {y), 

m=—l 

where we take x,y on the unit sphere S 2 , (x,y) is the standard scalar product on K 3 . 
Let us now define the kernel 



AM*>V))= E b 2 (^-)L l {{x,y)l 

[B3- 1 ]<i<[Bi+ 1 ] 



[•] denoting integer part. Finally, we introduce cubature points, that is, for each I we 
consider {77:77 G A"/}, a finite subset of § 2 , and positive real numbers A I( > (the cubature 
weights) indexed by the elements 77 of Xi, such that 



V./' 



eJC h [fdn= E Vfa)- (2) 

J vex, 



It can be proved that, if X\ is a maximal e-net (see the definition below) with e ~ _B~ J 
for some j, then it constitutes a set of cubature points (see [19]). This will be our choice 
in the sequel. We will write Zj = X 2 [B3+i], Zj = ■ ■ • i £jk, ■ ■ •} and Xjk = X v in order 
to stick with the usual wavelet notation. 

We are now ready to introduce needlets, which are given by 

^jk(x):= E \/%kb(^-jLi((x,^ k )), j = l,2,...,£ jk eZj. 

[B3- 1 ]<Z<[B3 + 1 ] ^ ' 

Needlets enjoy a number of very important properties that are proved elsewhere and that 
we summarize in the following two propositions. Write as usual 



we have the following proposition. 
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Proposition 1 ([19]). (a) (Reconstruction) The family {ipjri)j&i,rieZj * s a tight frame, 
hence 

f( x )=Yl H (f,ipjk)L^)ipjk(x). (3) 
(b) (Localization) For any positive integer M there exists a constant cm such that: 

^ (x)l -irT^f^F' (4) 

where d(x 7 y) = arccos((a;,2/}) denotes the usual geodesic distance on the sphere. 

Tight frames in some sense are very close to orthonormal bases, as shown by (3); see [12] 
for further references and discussion. The good localization properties highlighted in (4) 
are very useful when dealing with missing observations or with statistics evaluated on 
subsets of the sphere (see below). 

Let us now focus on zero-mean, mean square continuous and isotropic random fields. As 
is well known, these can be expressed by means of the following spectral decomposition, 
which holds in the mean square sense: 

oo 

l — l Ira 

where {a/ m }z )rn , m = 1, . . . ,1 is a triangular array of zero-mean, orthogonal, complex- 
valued random variables with variance E|a; m | 2 = C;, the angular power spectrum of the 
random field. For m < we have ai m = (— l) m a;_ m , whereas a;o is real. Our random 
necdlct coefficients are then easily seen to be given by 

Pjk ■= I T(x)i}) jk {x) dx = v/Ajfc Yl h \w\ Tl ^ jk ^ € jk £ Zr 

JS2 [Bi- 1 ]<l<[Bi + 1 ] ^ ' 

Remark that in the previous sum only a finite number of non-vanishing terms appear. 
Hence we obtain 



l>0 



and 

Ei b2 Q>/ Bj ) c i L i((tj,k,tj,k>)) _ b 2 (i/B^cM(^k,^ k ')) 



COTT(/3j,k, Pj,k') 



Ei Vil/BWMl) Zi PQLmC&l + 1)/(4tt) 
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In particular E[/3| fe ] is equal to 

A*£^(4W«(1)- (5) 



z>o 



=73 



The following assumption is a mild condition on the regularity of the angular power 
spectrum. 



Condition 1 . Let us assume that 3a > 2 and an integer M > 3 such that 

BiJ' 



Ci = l- a 9j[-^), /or B^ 1 < I < B j+1 



W = 0, 1,...,M, sup sup k (r) (u)l < C r < oo. 

j l/B<u<B 

Remark 2. A typical example that fulfills Condition 1 is 

C l = G(l)l- a , sup u r \G {r) (u)\ <C r <oo, r = 0,l,...,M. 

u>l/B 

The following proposition provides the basic correlation inequality for the analysis of the 
statistical procedures in the sequel. 

Proposition 3 ([4]). Under Condition 1 we have, {or each positive integer M , 



CmB' 2 ^ 
{l + BH{i hk ^. k ,)Y 



Remark 4- The previous inequality is the basic ingredient for many of the results to 
follow, together with the geometric analysis in the next section. Inequality (6) entails 
that needlet coefficients are asymptotically (in the high-resolution sense) uncorrelated at 
any fixed angular distance. This makes consistent inference procedures viable, even in 
the presence of observations on a single realization of a random field. 



3. The geometry of Voronoi cells on the sphere 

In this section, we establish some properties of Voronoi cells on the sphere. These prop- 
erties will be instrumental for the correlation inequalities we shall need in our main 
arguments. Let us first introduce some notation. 

We start by defining the standard (open and closed) balls in § 2 as 



B(a,a) = {x, d(a,x)<a}, B°(a,a) = {x, d{a, x) < a}. 
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Also, if A C S 2 we denote by |^4| the two-dimensional measure of A. Of course, for the 
full sphere |§ 2 1 = 4n. Now let e > and xi, ■ . . , xn £§ 2 such that 

Vi^j, d(xi,Xj)>s 

and the set {x\, . . . , xn} = S e is maximal for this property, that is, 



VxeS 2 , d(x,E £ )<e and Vi^j, B[ x i} | j f\B\Xj^ 



We call such a set 5 e a maximal e-net. 

Of course, if e = 7t, obviously JV = 1, and as soon as e < tt, N >2. For instance, if 
5 < e < 7t, we can take TV = 2 and £i and £2 are two opposite poles; however, N = 3 is 
possible taking three points on a geodesic circle at distance The next lemma is a 
simple but useful result. 

Lemma 5. Let {xi, . . . , xjv} = S £ be a maximal e-net. Then 



A A 

±<jv<4^. 

S2 



(7) 



Proof. Let us first recall that Va £ S , VO < 77 < 7t, 

\B(a,rj)\ = 2n sin 6» d# = 27t(l - cost?) = 4tt sin 2 (77/2) (~ to; 2 for 77 -> 0). 
Jo 

More precisely, if < 77 < 2a < n, as i sin a > ^ , 

4 



.,1 o / sin a , 

77 2 -<7t7? 2 <\B(a,ri) \<KV 

7t V a 



(8) 



As 



we have 



(J 5(^,e) = S 2 , 



47t=|§ 2 |< \B(xi, e)\ < Ntxe 2 



and, as the B(xi, |) are disjoint 



2 n 



<4n. 



□ 



Subsampling needlet on the sphere 445 
Remark 6. Actually, it is easy to see that 

1 <A< 1 



sin e/2 sin e/4 

so that 

11 e 

- < s— r- < Asm 2 - < 1. 

4 4cos^e/4 4 

We can now introduce, for all Xi 6H E , the Voronoi cells: 

Definition 7. Let S e be a maximal e-net. We define the associate family of Voronoi 
cells: 

S{xi) = {x £ S 2 ,Vj ^ i,d(x,Xi) < d(x,Xj)}. 

Clearly 

Now observe that \/a^=b and the set {x G S, d(x, a) — d(x, b)} is a geodesic circle on the 
sphere. Let us denote this circle by C a > The set: {x 6 S 2 ,d(x,a) < d{x,b)} = D(a,b) is 
the half-sphere defined by this geodesic circle containing a. We have 

S(xi) = P| D(xi,Xj), 

so S(xi) is the intersection of a finite number of half-spheres. Actually, 

S(xi)= P| D(xi,Xj) 

and, more precisely, 

S(xi)nS(xj)^0 d(xi,Xj) < 2e. 

The following lemma provides a bound on the number of points of a maximal e-net that 
lie within a distance 2e from a fixed center. 

Lemma 8. 

Card{j ^ i, d(xj,Xi) < 2e} < 6tt 2 . 
Proof. Let K(e,Xi) = {xj ^ x i} d{x 3 -,, Xi) < 2s}. Clearly, as 

K(e,Xi) C B(xi,2e)\B(xi,e), 
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we have, if Xj G K(e,Xi), B(xj,^) C B(xi, ^) \B(xi, §). Hence, since the balls B(xj, §) 
are disjoint, 



E 



S [ Xj ,- 



< 



B( ) \B[ x 



For < n < i] < 7t, we have 

f 

|-B(a,7?) \B(a,/i)| = 27t / sin d9 = 27t(cos — cos rf) 
= 47tsin 

Hence 



sin 



<^ 2 -M 2 )- 



7t 

V< 7T 



-(r/ 2 - M ) < |B(a, 77) \ B(a,/z)| < tt(t7 2 - /i ). 

7T 



Now, if r\ — > and < < 77, 



|B(a,7 7 )\B(a, / x)|~7t(7 7 2 -/i 2 ) 



So, by (9) and (8), 



which implies 



£ \ 2 4 



card^e,^))! - J -<tt( ( - 



card(^(e,Xi)) <6tt 2 , 



(9) 



(10) 
(11) 



as required. 



□ 



Of course 

|J S( Xi )= |J B(x l ,e) = E> 2 . 

It is obvious that if x, y € S(xi) then the portion of geodesic circle joining x to y is inside 
S{xi). 

Definition 9. Two Voronoi cells S{xi) and S(xj) with i ^ j are said to be adjacent if 

S(xi)nS(xj)^0. 

Remark 10. From Lemma 8 it is clear that there are at most 6n 2 adjacent cells to any 
given cell. 
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With the previous results at hand, we have now the ingredients to define our statistics 
on suitable subsets of S 2 , as explained in the following section. 



4. Asymptotics for needlet functionals 

In Baldi et al. [4] we discussed the asymptotics of several statistics of the form 



T< J 9] '^irE^^)' where 



Here Aj denotes the number of cubature points and H q as usual are the Hcrmitc poly- 
nomials of order q, (see, e.g., Surgailis [27]): 

i^) = (-l)V^£_e- 2 / 2j 

that is, Hi(u) = u, H2(u) = u 2 — 1, H^u) = u 3 — 3it, .... It was shown that ET^ = 
and, as j ^ oo 



r (9) 
3 



:iV(0,l), where £ f> :=E[(TV" ) 2 }. 



Indeed it is also possible to show stronger results; for instance, the joint convergence over 
different values of q. 

Our purpose here is to provide estimates for the normalizing variance . Our idea 
is to "subsample" our statistics by evaluating them on distinct Voronoi cells. By means 
of the asymptotic properties of needlet coefficients, we shall show that each cell will 
provide asymptotically an independent realization of the same limiting distribution as the 
statistic on the whole sphere. It will then be immediate to exploit these results to estimate 
by bootstrap our limiting distribution. More precisely, fix B > 1, let j,r be two non- 
negative integers such that B~ J < B~ r /A and write "B nB -j := {xj t i, . . . , Xj y A, }, ^nB- r '■= 
{av,i, . . . ,x rt A r } for corresponding sequences of maximal e-nets with e = nB -3 and = 
tcB~ t , respectively. Let us also define 

N a ;rj :=Caxd(S(x rta )r\E„ B -j), A r = Caxd{E nB -r}, 

where S(x r , a ) are the Voronoi cells of the 7tS _r -net E nB - r . In other words, N ra; j rep- 
resents the number of points in the B~ 3 -net that fall inside the Voronoi cell of E nB -r 
around the point x ra (note that x r ^ a need not belong to E nB -j). A r denotes the cardi- 
nality of such Voronoi cells. In particular, 5Zo=i N a -,rj = Ar 
We focus on the triangular array 

T atrj :=-^= J2 H «(M, a=l,2,...,A r , (12) 

v 7V ^- k£3nB _. nS(Xria) 
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for a fixed positive integer q. We define cr| := Vax{H q (/3jk)} , the variance of each sum- 
mand in (f2), and S a;r j := Var{r a;r j} the variance of the whole sum. In the sequel, we 
drop the index q whenever this causes no ambiguity. In particular, for r = the unique 
Voronoi cell is S 2 itself and we write Sj := £ a: oj- We prove the following: 

Theorem 11. Assume that the field is Gaussian and that Condition 1 holds. 
(a) As j — * oo 

lim — y = 1, in probability, 

j-*OQ Oj 



where 



8? 



(b) Forr,j such that 



Sj is consistently estimated by 



A 3 



- + '-0, (13) 



S i=i-E r ^> ( 14 ) 



A r 

a — 1 

that is, 

% 

lim — - = 1, in probability. 
(c) For r (and hence A r ) fixed and j — > oo, we have 

l;rj r Arjrj 1 



iV A ,.(0,/), asj^oo. (15) 



(d) As r — > oo, and £ — ► oo, we /lave 
r 

1 rn.: 



{^«}a=I,...,A r , Z a ±NID(0,l), (16) 

^rj J a=l,...,A T , 

1 -{£a}a=l,...,^, Z o ^JVJ£>(0,l), (17) 

=l,...,A r 

the convergence taking place in the space of sequences t°° endowed with the standard 
topology. 
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Remark 12. We note here that (a) and (b) of Theorem 11 provide some sort of boot- 
strap/subsampling approximation for the sample variance for the statistics 



1 

fees. 



thus eliminating the need for parametric assumptions on nuisance parameters or lengthy 
Monte Carlo simulations. In particular, for q — 1, point (a) yields a consistent estimate 
of the normalizing factor jj defined in (5). 

Condition (13) of Theorem 11 is clearly reminiscent of standard bandwidth assump- 
tions in nonparamctric estimation. It ensures that the number of pixels on which needlet 
coefficients can be evaluated grows more rapidly than the number of Voronoi cells or, 
in other words, the number of observations within each Voronoi cell goes to infinity. 
Points (c) and (d) entail that statistics evaluated on (possibly adjacent) Voronoi cells are 
asymptotically i.i.d. 

In the sequel, we take r to be an increasing function of j, but we refrain from writing 
Tj to simplify notation. 

The proof of Theorem 11 requires a very careful bound on cross-correlation of our 
statistics over different Voronoi cells. Such bounds rely very much on the following result, 
the proof of which is quite delicate and collected in the Appendix. 

Proposition 13. Let < e < 8/ A and let Sa be a maximal S-net and S e be a maximal 
e-net. For x a £ S^, let S(x a ) denote the Voronoi cells centered at x a for the S-net Hj. 
Let N a -s £ = Card(S'(x a ) n S e ) and assume that e = nB^ . Then, if x a , Xb S Sj, x a ^ Xb 
and M > 3, there exists C > verifying the following property: 



e ns(x a ) u£~. e ns(x b ) 

Proof of Theorem 11. From Proposition 13 it is immediate to get the following in- 
equality, which we shall exploit several times in our arguments below. For some C > 0, 
we have 

(18) 

<C{j ~r)B- {j - r) logB. 

(a) The proof of (a) is very similar to the proof of b), indeed slightly simpler, and 
hence omitted. 

(b) It is readily checked that lim^oo EjEj/Ej} = 1. The idea of our argument is then 
to establish 

lim sup Var(^-^ij = 0(1), (19) 
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lim sup max 



J— >oo 



Cov 



■ ra;j rb\j 



= 



(20) 



and then to use these results to conclude the proof by noting that 



lim sup Var{5l| = lim sup -^ V Covj-^-, 

j_oo Lij J j^co A r a b=1 L 



1 Ar rr 2 i / 1 Ar rr 2 r 2 

limsup-l^Var ^ +0 (± £ Cov ^ ^ 



To make this argument rigorous, we start by noting that 

rr 2 r 2 

Covl ^l^J^Lt 

1 

ra; 7'^V r&; 7 



(21) 



E 



E{# 9 (& fcl )H q (0 jk2 )H g (0 jk3 )H q (0 jk4 ) } 



fa,beH KB _j ns(i„) 



Let pj (fcx , fe) = con(j3jki 1 Pjk 2 ) = corr(/3jfe 1 , /3j-fc 2 ) . By the diagram formula (see, e.g., Sur- 
gailis [27]) and with standard manipulations it can be shown that (21) is bounded by 



' J J k 1 ,k 2 &^„ B -jr,S(x ra ) 



x pf (k u ki)pf {k 2 ,k 3 )pf {k 2 , k A )pf (fc 3 , fed) I 



(22) 



for some K > and for all choices of non- negative integers q% , . . . , q@ such that qi + ■ ■ ■ + 
qe = 2q and (q 2 V q 3 ), (q^ V q§) are strictly positive. Rearranging terms and recalling that 
|Pj(v)l 1) we obtain 



(22) < 



K 



' J J k 1 Me=. nB - ] nS{x ra ) 



\Pj(ki,k 3 )pj(k2,kA) 



N ra -jN r b-j 



E \Pj{ki,k3)\ 
fciea 7tB _ i ns(x ro ) 
^e=„ B -jnS(i r! ,) 
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< 



E 



ft 3eH„ B _3nS(a: r! ,) 

< C(j - r) 2 B~ 2 ^-^ log 2 5, 

in view of Proposition 13 and (18). Thus (20) is established. The proof of (19) is similar. 
It follows that 



lim Vari =^ 1 = lim Var< — 

1 Ar f 

= lim —77 >^ Var< 
i-oo ^ 1 
' a— 1 v 



ray 



r 2 ■ 



1 Ar f 

lim — Cov<^ 



V 2 V 2 

ra;j rb:J 
Ej Sj 



< o 



lim sup max 



, Ka,b<A, 



Cov 



r 2 r 2 

ra;j ro;.7 



:o(l), asj^oo, 



which completes the proof of (b). 

(c) Here we recall that A,- = A is fixed, that is, we are focusing on a finite number of 
subsets of the sphere. In these circumstances, the argument is very similar to the proof 
of Theorem 9 in Baldi et al. [4]. We use the Cramer- Wold device and thus focus on the 
linear combination 



-= W 3^ja N I 0, W ja J ■ 

V U i o=l V o=l / 



where Wj' a £ R for all j = 1, 2, . . . , a = 1, . . . , Aj. It is obvious that 



On the other hand 



Var< 



£a=l W Ia E i 0=1 



Ea=l W %^j a^b 



^2wj a w jb E{Tj a T jb } -> 1, 
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as j — > oo, in view of (20). Also, let us denote by cum p (X) the cumulant of order p 
of the random variable X; to establish the central limit theorem, we resort to a recent 
result by Nualart and Peccati [21], where it is proved that convergence to zero of the 
fourth-order cumulant is a sufficient condition for asymptotic Gaussianity in the Gaussian 
subordinated case we are considering here (see also DcJong [7] for a similar approach 
with multilinear forms in i.i.d. sequences). Again by the diagram formula for cumulants, 
we obtain easily for p — 4, 



3 a 3 6 J c J ^ 



< 



a,6,c,d=i W^i v 2 ^' y^i V^3 

KA$ 

\J N ra -jN r b;jN rc -jN r d-j 

H q (fi jkl ) H q (P jk2 ) H q (p jk3 ) H q (f3 jk J 



y^j y^o y^j 



x cum 

k 4 ee- 7xB - j nS(x rd ) 

y N rai jN r b-jN rc -jN r d-j 

x fc 2 )| |/0j- (fc 2 , fc 3 )| 1/0^(^3, A; 4 ) | (fc 4 , /ci) 

fc46H 7lB _ i nS(2: rti ) 



fc2e3„ B - 3 nS(x rb ) 



1 V^W^-,,^^,^, (1 + B^d{^ jM )) 
k i e F z nB - j nS(x Td ) 

= o(l), asj^oo, 

again in view of Proposition 13 and (18). The central limit theorem is then established. 

(d) Fix a finite subset D c N of cardinality A and label its elements a = 1, . . . , A. It 
is obvious that as j — ► oo 

r rfcw -\ _/r™a . d ^iiv(o,i), 



J t£D I V J o=l,...,i4 
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by exactly the same argument as in (c). Because D is arbitrary we have thus established 
convergence in the finite-dimensional distributions sense of the sequence {r ra;j -/->/S^} ae i!j 
to the sequence {Z a } aS N- It is a standard result that in R°° the finite-dimensional sets 
are a determining class (see Billingsley [5], page 19), so that weak convergence is actu- 
ally equivalent to convergence of the finite-dimensional distributions. The argument for 

{r r a;j / y^j}aen is entirely analogous by means of Slutzky's lemma. 

□ 



5. Discussion and statistical applications 

The results in the previous section lend themselves to several applications for the sta- 
tistical analysis of spherical random fields, in particular for CMB data. Here we briefly 
discuss some examples; we do not provide a complete discussion or application to real 
data, as we prefer to defer this more detailed investigation to future works. Refer also to 
[9, 22] for numerical evidence and applications of needlets to CMB data. 

An immediate application of the results in Theorem 11 concerns studentization. Indeed, 
in Baldi et al. [4], several statistics were proposed, based on needlet functionals, to test 
for goodness of fit, Gaussianity and isotropy of CMB data. We remark that these three 
issues are very widely studied in the huge satellite experiment collaborations on CMB, 
such as WMAP and Planck. These statistics could all be expressed as linear combinations 
of {r ra; j} in (12), and asymptotic Gaussianity was established under the null, provided 
the normalizing variance could be taken to be known. The results of the previous section 
can then be immediately exploited to provide estimates of the limiting variances, thus 
making feasible testing procedures with a standard asymptotic distribution. In particular, 
Theorem 11, part (d) shows how nonlinear statistics can be studentized by the estimates 
Sj given in (14) without affecting the limiting distribution. 

Other possible applications of our results relate to testing for isotropy of spherical 
random fields data. Indeed, in cosmological applications an issue that has drawn an 
enormous amount of attention over the last three years is the possible existence of sta- 
tistical asymmetries in the behaviour of CMB radiation. The assumption of an isotropic 
universe is very much embedded at the roots of cosmology, in view of the so-called Ein- 
stein cosmological principle that the universe should "look the same" to any observer. 
However, quite unexpectedly, some evidence of statistical anisotropy is indeed present in 
the first releases of WMAP data (2003, 2006); see, for instance, Hansen et al. [10] and 
the references therein. Statistical procedures considered so far have led to inconclusive 
results. This mixed evidence has sparked an enormous amount of further empirical re- 
search, as such asymmetries may entail profound consequences in fundamental physics 
(as mentioned in the Introduction) . It is therefore of great importance to devise new sta- 
tistical tests that can exploit as efficiently as possible the available data. In view of their 
double localization in real and harmonic space, needlets emerge as natural candidates 
to build such procedures. We stress, in fact, that the knowledge of the scales where the 
asymmetries might lie would provide essential information toward their understanding. 
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We mention here that many other new wavelet-related systems have recently been intro- 
duced to deal with anisotropic features, such as curvelets and other forms of directional 
wavelets (see Jin et al. [13], Starck et al. [26], Vielva et al. [30] and Wiaux et al. [32]). 
We view these important approaches as complementary to ours and we leave for future 
research the investigation of interactions between these lines of research. 
We suggest the following procedures: For fixed q and A, consider 

r 



Sj := sup 

a=l,...,A, 

under the null of Gaussianity and isotropy, we have from Theorem 11, part (c): 

lim Px{Sj < x} = {$(x) - $(-.t)} a '\ 

3— >°o 

where $ denotes the standard cumulative distribution function of a Gaussian random 
variable. Under the alternative, wc would expect to obtain unusually large values over 
some regions of the sky where data are generated according to a different model, and thus 
we expect the procedure to have good power properties against suitable alternatives. 

An alternative approach can be envisaged as follows: In many applications, a possible 
direction for the asymmetries is easily conjectured. For instance, with CMB data, natural 
frames of reference are provided by either the ecliptic plane (i.e., the subspace approx- 
imately spanned by the planetary orbits around the sun) or the galactic plane, that is, 
the plane approximately defined by the location of the Milky Way. Asymmetries in these 
directions would lead presumably to explanations that do not directly involve the CMB 
itself, but rather other astrophysical entities of a non-cosmological nature. These forms 
of north-south asymmetries can be readily tested by means of the procedure that we 
describe below. Let 

H q (J3 jk ), a = 1,2, 



where xi,X2 £ S 2 denote the north and south poles in the suitable frame of reference 
(e.g., the so-called galactic or ecliptic poles). From Theorem 11, part (c), we have that 

Tj := - Tj 2 . 3 } 2 ->d Xi under H , 

thus immediately providing threshold values with an asymptotic justification. 

More sophisticated approaches are indeed possible. One idea is to exploit the asymp- 
totically i.i.d. behaviour of ncedlet coefficients to transform these statistics into an ap- 
proximate sample of spherical directional data, and then use the rich machinery of testing 
for uniformity that has been developed in these circumstances (see Mardia and Jupp [14], 
Pycke [24] and the references therein). A possible idea is to proceed with a hard thresh- 
olding of needlet coefficients and then focus on the directions that correspond to selected 




Subsampling needlet on the sphere 



455 



coefficients. More precisely, we can define the new data set 

{x rj , T eS 2 :|/? jfe | >Tj}, 

for a suitable choice of the thresholding sequence {tj}. Under the null, in view of the 
limiting properties of {Pjk}, the sequence {x r j [T } is approximately uniformly distributed 
on the sphere, a conclusion that can easily tested by a variety of well-established pro- 
cedures. A rigorous investigation of this proposal, however, is beyond the scope of the 
present paper and will be deferred to future research. Likewise, applications of the ideas 
in this section to CMB data from WMAP are currently in progress and will be reported 
elsewhere. 



6. Appendix: Proof of Proposition 13 

For notational simplicity and without loss of generality, throughout this Appendix we 
take B = 2. The proof requires several steps. The first part, which is relatively simple, 
refers to circumstances where the distance between Voronoi cells is larger than their 
radius. 



6.1. First case d(S(x a ), S(x b )) > 6 

Under this hypothesis, and due to the next lemma, we have 
i 



E E 



y/N a ,ssN b , Se c _^. . c Jrt. (l + 2J - d(u,t;)) 



i / r\ 2 / \ M / \ M-2 



Lemma 14. Let < e < | and let £5 be a maximal 6 -net and S e bea maximal e-net. 
Let x a G Si and let S(x a ) be the corresponding Voronoi cell. Then 

1 „ . _ , _ o/5 N 2 



-J ^<Card(5(x a )nS £ )<27t 2 ( - ) . (23) 



Proof. As B(x a , |) C S(x a ) CB(x a ,6), we have 



|J B[u,-)cB(x a ,S+- 

«eS(i„)nH e 



and 



B(x a , 5 --e)(l |J B(u,i 
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In view of (8), it follows easily that 

2 / \ 2 
Card({u 6 S(x a ) nS e })^- <tt((5+| 



Card({it e S{x a ) n S e )}) < 2tt 2 



^ 2 



Hence 



and moreover 



hence 

/ 2 i ^<^^-l) <Card({uG5(x a )nS £ }). □ 
6.2. Second case d(S(x a ), S(x b )) < 6 

Here we face the situation where we have neighbouring Voronoi cells; the corresponding 
covariances are clearly harder to bound and we shall first introduce several lemmas. 

Lemma 15. Let e ~ 2~ J and M > 3. Let x a ^ Xb and x ai x b £ H{. Let u G S e n S(x a ) 
be fixed. Then 



-(tj- 2 ) <Card({uGS*(a; Q )nS e })7t£ 2 , 



A 2 1 4/5 



<c M - 



(l + 2id(u,S(x b ))) M - 2 - 
Cm can be chosen equal to 27t3 2A/ ~ 1 . 

Proof. By straightforward manipulations, we obtain the bound 



b) (l + 2^))M 

= y I / I dx 

— V f — 
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by (8) and as by the triangle inequality d(u,x) < 2d(u,v) for x <S B(v,e/2). Clearly 

f 2 M 

- J s/2 (1 + V6) M ~ J e/2 (V6) M ~ M-2 
which implies 

A{u)<tx2 2M -\ 
s 

The previous inequality for d(u, S(x a )) < — yields immediately 

3 2M-1 1 



(l + Vd(u,S(x b )))M-f 

On the other hand, if d(u,S(x a )) > 



2 

f 2 M 

f 2 M 
< / ([x 

Jd{u,S(x a )) (1 + 2 J t>) Jd(u,S(x a )) 

<2n2->j^(Vd(u,S(x a ))-( M -V <2n2-\yd(u,S(x a ))))- {M - 2 \ 



whence we get 



A{u) < 2n{Vd{ Ul S{x b ))- {M - 2) < 27t3 m " 2 - 



(l + 2Jd(u,S(x b ))) M - 2 ' 
We are thus able to conclude that 



V 1 <2 



2M-1 



b) (l + 2id(u,t0)" " (l + 2^(u,5(x fc )))M- 
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Using the previous Lemma 15, we have: 



drf 1 y y 

. /AT - AT. - ^ ^ 



' »£= E nS(i„)i)£E e nS(ib) 



< 1 2 7 t3 2M - 1 y 



V^a^fe ue5 ^ s(iJ (1 + 2*d(u, 5(.T fc )))^-2 ' 

Now as d(S(x a ), S(xb)) < 5, let us observe that if w £ S e n 5(x Q ) 

d(u, S(x 6 )) < d(S(x a ), S(x b )) + 26< 3d. 
Some computations yield, using Abel's formula of summation by parts: 

27t3 2M ~ 1 Card{M G S £ n S(x a ),le < d(u, Sfo)) < (I + l)e} 

2tt3 2A/ - 1 Card{ueS e nS(x a ),(Z-l)e<d(u,S(x b ))<Ze} 
2tt3 2A/ - 1 f / e x M " 2 



< 7= AT =jf Cardie H e nS(a; ) ) d( U) 5(a; 6 ))<3* + e} 



Z (]ia _ (1 + ^/- 2 ) Card{w G ^ £ n S(x a ), d(u, 5(x 6 )) < Ze} J 



l<Z<3<5/e 



(57 I V 35 



\ I _ 2 

zJ ;a/-i Card { w G 2 e n Sfca), d(u, 5(x b )) < Is} 

l<l<36/e 



/ \ M-2 

i 



<167t 5 3 A/+1 (J 

+ 87T 3 3 2A/-i/£\ y M^Caxd{ueE £ nS(x a ),d(u,S(x b ))<le}. 

^ ' l<;<35/e 

Let us introduce a further auxiliary result. 
Lemma 16. For d(S(x a ), S(x b )) < S, e < 5 /A, 

Card{u G S e n S(x a ), d{u, S(x b )) < Ze} <d(l + l) 2 + C* 2 -(Z + 1) 

e 



(with, e.g., C\ = 6n 4 and C 2 = 547t 4 / ). 
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Proof. Using Lemma 8, if dA denotes the boundary of the set A, we have 

A I 



d(s( Xb )) = U r m , m < 6ti 2 



m— 1 

where r m = [a m , Om+i] is a portion of a geodesic circle T m C C Cm := d(B(c m , 5) of length 
less than diam(5(x&)) < 2(5, c m denoting its center. Let u £ S e n S(x a ), d(u, S(xb)) < Is- 
Certainly 

d(u, S(xb)) = d(u,w), w€d(S(xb))- 

We split the proof into two parts according to whether w is a corner point or not. 
• If w is a corner point a m , then 

it£S e , d(u, a m ) < le < 36 

and 

Card{u G 3 e , d(it, a m ) < le} < \n 2 (2l + 1). 
In effect, for such u we have 



So if k m = Cardjw <E S e , d(u,a m ) < le}, as / > 1, using (8) 

fc 4(l) 2 ^0 + ^) 2£2 ' 

hence 

fc m < i 7t 2 (2Z + l) 2 < 7T 2 (Z + l) 2 . 

• On the other hand, let us focus on w G]a m ,a TO +i[, 

d(u, S(xb)) = d(u, w) = d(u, C Cm ) < 36. 

So we have that 

and also to the portion of the half sphere B{c mi n/2) , which is between the two geodesies 
joining c m to a rn and c rn to a m +i- Let us call this set G(c m ;a m ,a m+ i). 
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Therefore {u e S e n S(x a ), d(u,S(xb)) < le} is contained in the union of the M sets 
{« £ 2 E fl B(a m ,Ze)} and in the union of the M sets {11 6 H £ fl G(c m ;a m ,a m+ i) n 
B(c m , ^) \ B(c m , ^ — Is)}. Let us evaluate the cardinality of each of these sets. 

Let us consider the set {u £ E e n G(c m ;a m , a m+ i) n B{c m , f )\-B(c m , § — ie)}- It is a 
simple observation that for a point m of this set, 

^(^|)^(c ro) f + |)y(c m ,f-(/ + i) £ ) (24) 

(recall that the -B(it, |) are disjoint), and also 

where [a m , a m +i] is a portion of geodesic circle (L rn C C Cm containing [a m , a rn +i] of length 
28 + 2a, with a < <5. Indeed, such a point w is on a geodesic joining c TO to w £ ]a m , a r „+i [. 

Let us recall that d(c m ,w) = 5 and thus as d(u,w) < 3(5 and if 5 < j|, certainly, the 
ball B(u,e) does not contain the point c m (as e < |). 

Hence, let us consider the geodesic tangent [c m ,w], where v belongs to the boundary of 
the ball. Let us also consider the spherical triangle (c m ,u,v). The a that we are trying 
to bound is the angle at the vertex c m . 

By a standard spherical trigonometric formula, we have: 

sin a sin v 



sinc?(u,u) sind(u,c m )' 

where v denotes the angle at the vertex v in the spherical triangle (c m , u, v). This formula 
translates here as: 



sin a sin7r/2 



Hence, as e < | , 



sine sin(7t/2 — d(u, w)) cosd(u,w) 



sine sm5/4 
sin a < < — < smo. 

cos 3o cos 3<5 



The last inequality is easy to check for 8 < j| . 

We are now in the position to conclude the proof of Lemma 16. 

Let us first take 5 < jk. So, if k m = Card{u £S £ fl G(c m ;a m ,a m+ i) n B(c m ,^) \ 
B(c m ,% -le)}. By (24) 



7l\2 



Cm, ^ + |)\- B ( Cm ' \ _ y + \ )' 
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So 



4 e 



n\2 



AS 

< n ^1 

~ 2tt 



7T £ 

2 + 2 



71 



2rt 



7T(Z + l) £ + -e 2 - Z + - e 



1 



< 2tt5(Z + l)e. 



So 



fc m < -2tT(/ + 1). 

On the other hand, let us suppose instead that j| < 6 < n. The computation of the 
cardinality K = Card{u 6S £ n B(c m , f ) \ B(c m , f — le)} is much simpler: By (24) 



K- 



7T\2 



< 



B[c m , 



7T £ 

2 + 2 



B c 



71 



which leads to 



4(1' 



7t e 

2 + 2 



<4 



<7T 2 (Z + l)e 



and therefore 

K < -n 3 (l + 1) = — 127t 2 (Z + 1) < -127t 2 (/ + 1). 

£ 12 £ £ 

By the previous computations and recalling that the number of adjacent cells is < 67t 2 , 
we have completed the proof of Lemma 16. □ 



Now we can conclude the proof of Proposition 13 by noting that 

, M-2 

W < 16tt 5 3 m+1 f £ - 



87T 3 3 M " 2 



M — 2 

^2 im-i Caxd{uGE e nS(xi), d(u,S(x b )) < le} 



< 167T 



A/-2 



l<l<3S/e 



87T 3 3 M-2 



l M 

£ 

6 



E ^<Mi + D a 

l<Z<3<5/£ 



l<Z<3<5/e 
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and if M > 3, 



as well as 



l<i<3<5/e v 7 l<Z<3<5/e 



l<;<35/e l<i<35/e 



References 

[1] Antoine, J.-P. and Vandergheynst, P. (1999). Wavelets on the 2-sphere: A group-theoretical 
approach. Appl. Comput. Harmon. Anal. 7 262-291. MR1721807 

[2] Antoine, J.-P., Demanet, L., Jacques, L. and Vandergheynst, P. (2002). Wavelets on the 
sphere: Implementation and approximations. Appl. Comput. Harmon. Anal. 13 177- 
200. MR1942741 

[3] Baldi, P., Kerkyacharian, G., Marinucci, D. and Picard, D. (2008). High frequency asymp- 

totics for wavelet-based tests for Gaussianity and Isotropy on the Torus. J. Multivariate 

Anal. 99 606-636. MR2406074 
[4] Baldi, P., Kerkyacharian, G., Marinucci, D. and Picard, D. (2006). Asymptotics for spherical 

needlets. Ann. Statist. 37 1150-1171. 
[5] Billingsley, P. (1968). Convergence of Probability Measures. New York: Wiley. MR0233396 
[6] Cabella, P., Hansen, F.K., Marinucci, D., Pagano, D. and Vittorio, N. (2004). Search for 

non- Gaussianity in pixel, harmonic and wavelets space: Compared and combined. Phys. 

Rev. D 69. 

[7] DeJong, P. (1990). A central limit theorem for generalized multilinear forms. J. Multivariate 

Anal. 34 275-289. MR1073110 
[8] Dodelson, S. (2003). Modern Cosmology. London: Academic Press. 

[9] Guilloux, F., Fay, G. and Cardoso, J.-F. (2007). Practical wavelet design on the sphere. 

Appl. Comput. Harmon. Anal. To appear. Available at arXiv:0706.2598vl. 
[10] Hansen, F.K., Banday, A. J., Eriksen, H.K., Gorski, K.M. and Lilje, P.B. (2006). Foreground 

subtraction of cosmic microwave background maps using WI-FIT (wavelet based high 

resolution fitting of internal templates). Astrophysical J. 648 784-796. 
[11] Hansen, F.K., Cabella, P., Marinucci, D. and Vittorio, N. (2004). Asymmetries in the local 

curvature of WMAP data. Astrophys. J. 607 L67-L70. 
[12] Hernandez, E. and Weiss, G. (1996). A First Course on Wavelets Florida: CRC Press. 

MR1408902 

[13] Jin, J., Starck, J.-L., Donoho, D.L., Aghanim, N. and Forni, O. (2005). Cosmological non- 
Gaussian signature detection: Comparing performance of different statistical tests. 
EURASIP J. Appl. Signal Process. 15 2470-2485. MR2210857 
[14] Mardia, K.V. and Jupp, P.E. (2000). Directional Statistics. Chichester: Wiley. MR1828667 
[15] Marinucci, D. (2006). High resolution asymptotics for the angular bispectrum of spherical 
random fields. Ann. Statist. 34 1-41. MR2275233 



Subsampling needlet on the sphere 



463 



[16] McEwen, J.D., Vielva, P., Hobson, M.P., Martinez-Gonzalez, E. and Lasenby, A.N. (2007). 

Detection of the integrated Sachs-Wolfe effect and corresponding dark energy con- 
straints made with directional spherical wavelets. Monthly Notices Roy. Astronom. 
Soc. 376 1211-1226. 

[17] McEwen, J.D., Hobson, M.P., Lasenby, A.N. and Mortlock, D.J. (2006). A high-significance 

detection of non-Gaussianity in the WMAP 3-year data using directional spherical 

wavelets. Monthly Notices Roy. Astronom. Soc. 371 L50-L54. 
[18] Moudden, Y., Cardoso, J.-F., Starck, J.-L. and Delabrouille, J. (2005). Blind component 

separation in wavelet space. Application to CMB analysis. EURASIP J. Appl. Signal 

Process. 15 2437-2454. MR2210855 
[19] Narcowich, F.J., Petrushev, P. and Ward, J.D. (2006). Localized tight frames on the sphere. 

SI AM J. Math. Anal. Appl. 38 574-594. MR2237162 
[20] Narcowich, F.J., Petrushev, P. and Ward, J.D. (2006). Decomposition of Besov and Triebel- 

Lizorkin spaces on the sphere. J. Fund. Anal. 238 530-564. MR2253732 
[21] Nualart, D. and Peccati, G. (2005). Central limit theorems for sequences of multiple stochas- 
tic integrals. Ann. Probab. 33 177-193. MR2118863 
[22] Pietrobon, D., Balbi, A. and Marinucci, D. (2006) Integrated Sachs-Wolfe effect from the 

cross correlation of WMAP 3 year and the NRAO VLA sky survey data: New results 

and constraints on dark energy. Phys. Rev. D 74 043524. 
[23] Politis, D.N., Romano, J. P. and Wolf, M. (1999). Subsampling. New York: Springer. 

MR1707286 

[24] Pycke, J.-R. (2007). A decomposition for invariant tests of uniformity on the sphere. Proc. 
Amer. Math. Soc. 135 2983-2993. (electronic). MR2317977 

[25] Sanz, J.L., Herranz, D., Lopez-Caniego, M. and Argueso, F. (2006). Wavelets on the sphere. 

Application to the detection problem. In Proceedings of the 14th European Signal Pro- 
cessing Conference (EUSIPCO 2006) (F. Gini and E.E. Kuruoglu, eds.). 

[26] Starck, J.-L., Moudden, Y., Abrial, P. and Nguyen, M. (2006). Wavelets, ridgelets and 
curvclets on the sphere. Astronom Astrophys 446 1191-1204. 

[27] Surgailis, D. (2003). CLTs for polynomials of linear sequences: diagram formula with il- 
lustrations. In Theory and Applications of Long-range Dependence 111-127. Boston: 
Birkhauser. MR1956046 

[28] Varshalovich, D.A., Moskalev, A.N. and Khersonskii, V.K. (1988). Quantum Theory of 
Angular Momentum. New Jersey: World Scientific. MR1022665 

[29] Vielva, P., Martinez-Gonzalez, E., Barreiro, B., Sanz, J. and Cayon, L. (2004). Detection of 
non-Gaussianity in the WMAP first year data using spherical wavelets. Astrophysical 
J. 609 22-34. 

[30] Vielva, P., Wiaux, I., Martinez-Gonzalez, E. and Vandergheynst, P. (2007). Alignment and 
signed-intensity anomalies in WMAP data. Mon. Not. R. Astron. Soc. 381 932-942. 

[31] Wiaux, I., Jacques, L. and Vandergheynst, P. (2005). Correspondence principle between 
spherical and Euclidean wavelets. Astrophys. J. 632 15-28. 

[32] Wiaux, I., McEwen, J.D., Vandergheynst, P. and Blanc, O. (2008). Exact reconstruction 
with directional wavelets on the sphere. Mon. Not. R. Astron. Soc. 388 770-788. 

Received June 2007 and revised August 2008 



